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Abstract 

We investigate spin-1/2 electrons with local Hubbard interaction and vari- 
able range hopping amplitudes which decay like sinh(K)/ sinh(Kr). Assum- 
ing integrability the Asymptotic Bethe Ansatz approach allows us to de- 
rive the generalized Lieb-Wu integral equations from the two-particle phase 
shift. Due to the nesting property there is a metal-to-insulator transition at 
Uc{k > 0) = 0^. The charge gap in the singular limit k = opens when the 

interaction strength equals the bandwidth, Uc{k = 0) = W > 0. 
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Typeset using REVT^ 



Exact solutions play a crucial role in our understanding of strongly correlated systems. 
The algebraic and the coordinate Bethe Ansatz constitute the main line of approach to inte- 
grable systems [jl|. Nevertheless, there are exactly solvable Hamiltonians whose eigenstates 
cannot be cast into the standard Bethe Ansatz form. A particular example is the popu- 
lar Calogero-Sutherland-Moser Hamiltonian |^,^ that has attracted recently a great deal 
of attention in connection with the universal properties of disordered systems and random 
matrix theory 0,^. 

Recently, Ruckenstein and one of us (F.G.) introduced a 1/r-Hubbard model that in- 
corporates long-range hopping and on-site repulsion, and that includes the lattice version 
of the Calogero-Sutherland-Moser model in the strong-coupling limit at half band- filling . 
The integrability as well as the structure of the eigenfunctions of this model are still poorly 
understood. The analysis of the two-body problem in the 1/r-Hubbard model indicates 
that the eigenstates are neither of the Bethe Ansatz form as in the Hubbard model, nor of 
the Jastrow-type as in the Calogero-Sutherland model Since the wave functions do not 
become plane waves even at large distances, Sutherland's Asymptotic Bethe Ansatz |]T],|^ 
cannot be directly applied to the 1/r-Hubbard model. 

In this article we introduce a model with variable range hopping that interpolates between 
the standard and the 1/r-Hubbard model. For half- filled bands and in the strong coupling 
limit it reduces to the antiferromagnetic 1/ sinh^(Kr)-Heisenberg or Inozemtsev model [§. 
This model is a special case of an exchange interaction model most recently explicitly 
solved in Ref. |jlO|]. The 1/ sinh(Kr)-Hubbard model is a straightforward but non-trivial 



generalization of the Inozemtsev model to an itinerant electron system [p!l|] . 



We consider a one- dimensional Hubbard model |jT2| for N spin-1/2 electrons with hopping 
amplitudes t{l — m) and on-site repulsion U, 

H = I] t(/ - m)ci+^c„^^ + UY1 , (1) 

ly^m.a I 

where the lattice sums on the ring run from —L/2 to L/2 — 1. The Hubbard model (|T]) is 
known to be exactly solvable in two cases: (i) t{l — m) = —t6i-m,±i for the standard Hubbard 



model [|T3| and (ii) t{l — m) = — 2t(— 1)' "^[d{l — m)] ^ for the 1/r-Hubbard model where 



d{l — m) = (L/tt) sin[7r(/ — m)/L] denotes the chord distance between the sites / and m. 
Henceforth we will use t as our energy unit. Here we introduce a model on the infinite chain 
{L —>■ oo) with hopping amplitudes t{l — m) = — i sinh(K)(— I)'"™'/ sinh[fi;(/ — m)] where 
controls the effective range of the hopping. 

The dispersion relation e{k) of our model (|l|) is given by e{k) = (— 2i) X^^i ^(^) sin(A;n). 
It is odd under parity, and can be expressed in terms of a logarithmic derivative of theta- 
functions. For n = one has e{k) = k with a discontinuity at the Brillouin zone boundary, 
while for k — >■ oo one finds e{k) = 2sin/c. e{k) is continuous for all k > with zeros at 
k = IT, k = 0. The maximum (minimum) at k = tt/2 {k = — 7r/2) for /t = oo is gradually 
shifted as function of k, to higher (lower) momenta until it reaches k = it {k = — vr) for k, = 0. 



The corresponding low-energy {g-o\ogj [T^) Hamiltonian involves left- and right-movers 



with different velocities. The limit k ^ is singular as the dispersion relation involves sud- 
denly only right-movers, while the left-movers' velocity diverges like 1/k, i.e., a discontinuity 
occurs in e{k) at the Brillouin zone boundary. Hence the physics of the metal-to-insulator 
transition for k, > and for k, = (1/r-Hubbard model) is completely different. At half- 
filling the distance = vrn between the two Fermi points becomes half a reciprocal lattice 
vector ("perfect nesting"). The bands for right- and left-movers effectively cross each other, 
and their degeneracy is lifted for all U > 0. The critical value for the transition is thus 
Uc{k, > 0) = 0^. In the chiral 1/r-Hubbard model where left-movers are absent one obtains 
a finite critical value, Uc{n = 0) = 2tt, see Ref. 0. 

We discuss the symmetric orbital part ?/'(xi,X2) of the two-particle wave function in the 
limit of a large system with open boundary conditions. In first quantization the Schrodinger 
equation for the wave function ^'(2^15 ^2) reads for xi X2 

E^{xi, X2) = F{xi, X2) + F{x2, xi) + t{xi - X2) {i^{x2, X2) - tpixi, xi)) , (2) 

while for xi = X2 we find {E — U)ip{xi,xi) = 2F(xi,xi). Here we defined F{xi,X2) = 
J2xj^x^,x2tixi - x)^p{x,X2). 
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We seek symmetric, spin singlet scattering solutions ip{xi <C X2) = e^i^^^'^+^'i^'i) — 
^ie{ki,k2)Qi{k2Xi+kiX2) ^ where 9{ki,k2) is the phase shift, and ki, ^2 are the quasi-momenta. 
Furthermore, we should choose ^/'(xi,X2) such that we can employ the trigonometric iden- 
tity sinh(2i) sinh(z2) = sinh(2;i — Z2)/ [coth(2;2) — coth(2;i)]. This naturally leads to 

^Xl,X2 



2 sinh[K(x2 — a^i)] 

_ ^-k{x2-xi) f j^^i{kiX2+k2Xi) _ j^^i{kixi+k2X2) 



k{x2-xi) ^j^^i(kiXi+k2X2) _ j^^i{k2Xi+k-i_X2)^ 

V5,„.,Ae*(^^+'=^)^^ (3) 



as choice for the wave function {B/A = e^^^'^'^''^'^)). It has precisely the form of Inozemtsev's 
two-magnon state 

The calculation of F{xi, X2) has to be done with care because of the infinite lattice sums. 
For xi ^ X2 the Schrodinger equation gives E = e{ki) + e(/c2), and 

{A + B){e{ki)-e{k2)) 



X = iA-B) + 



(4) 



2i sinh(fi;) 

while the equation for xi = X2 becomes {A — B)E — i{A + B){rj{ki) — rj{k2)) = X{E — U) 
with the abbreviation rjik) = 2 sinh(K) X]J^Li(~l)" cos(A;n) cosh(Kn)/ sinh^(Kn) = ri{—k). 
The phase shift is then found as 

'H{k2) - H{k^ 



e{kuk2 



-2tan"^ 



U/2 



(5) 



with H{k) = [-7]{k) + e{k){e{k) - U)/{2 sinh(/s:))] /2. 

Sutherland [|l]J^ observed that the two-particle phase shift is sufficient to obtain the 
(scattering) spectrum of a model, even if the wave functions cannot explicitly be constructed 
as long as the model is integrable. It is far from clear that the model (|l]) is integrable. 
Nonetheless we investigate the consequences of this conjecture. 

We employ periodic boundary conditions to quantize the pseudo momenta k^. The 
error introduced here is exponentially small, of order exp(— kL). It is clear how to set 
up the Asymptotic Bethe Ansatz equations for N — M up-spins and M down-spins in a 
straightforward generalization of the Lieb-Wu equations because H{k) plays the role 
of sin k: 



M 

Lkj = 2nlj + ^ e (2if(%) - 2Ap) (6a) 

(5=1 

N M 

-J2e (2A, - 2Hik,)) = 2n.J^ - ^ 9 (A„ - A;,) (6b) 

j=l 13=1 

for j = 1, . . . , N; a = 1, . . . , M, where Q{x) = — 2tan^-'^ {2x/U). Ij are integers (half-odd 
integers) for M even (odd), Ja are integers (half-odd integers) for N — M even (odd). 

We are interested in the ground state energy per site, E/ of the model in the thermo- 
dynamic limit. In this case the equations (^ can be transformed into integral equations for 
the densities p{k) = 1/ [L{kj^i — kj)], cr(A) = 1/ [(Aq,+i — Aq)] in the usual way: 

271 p{k) = l+H'{k) / d\a{\)K [H{k) - A; U/A] (7a) 

27ra(A) = f^^ dkp{k)K [H{k) - A; U/A] - f^' dX'a{X')K [A - A'; U/2] (7b) 

with K [x;y] = 2y/{x'^ + y^), and H'{k) = dH{k)/dk. The ground state energy density is 
calculated as E/L = Jq^ dkp{k)t{k). 

The integration limits Qi, Q21 Bi, and B2 are no longer symmetric around zero. They 
are determined by the particle numbers, Jq^ dkp{k) = N/L, j^^ dXa{\) = M/L, and the 
condition that a charge (spin) excitation at Qi {Bi) has the same energy as the corresponding 
excitation at Q2 (-82)- The latter conditions are most conveniently described in the pseudo 
particle picture ||15[. The pseudo particle dispersions for charge and spin for zero external 
magnetic field follow from 

edk) = e{k) + dXesiX)K [H{k) - A; U/A] (8a) 
27re,(A) = (^^ dkH'{k)t^{k)K [H{k) - A; U/A] - dX'ts{X')K [X - A'; U/2] . (8b) 

jQi JBi 

The integration bounds have to fulfill edQi) = ec{Q2), es(-Bi) = ^s{B2). The ground state 
energy and chemical potential can then be expressed in terms of the pseudo particle energies 
as E/L = J^^ dkec{k)/{27i), and p = MaxQ^<k<Q2 (^c{k). 

For zero external magnetic field, M = L/2, one finds B2 = 00 = —Bi. For half-filling, 
N = L, one further obtains Q2 = tt = — Qi- Therefore we can solve the integral equations 
analytically using Fourier transformation. With the help of the four functions 



,Im> 



dk 

-IT 2-K 

^ dk 

-n2^ 



,iu)H{k) 



e'''^^^'>e{k)H'{k) 



one may express the densities and the ground state energy as 

duo 



a{X) 



2'Kcosh{ujU/4:) 
1 



p{k) = - + H\k) 

2tt Jo 7t[1 + e'^^/^j 



cos{ujX)Jq{uj) + sm{ujX)jQ{uj 
doj 



cos{ujH{k))J^{uj) + sm{ujH{k))J^{uj) 



(9a) 
(9b) 



E/L 



°° 2duj 



(9c) 

-2Ji{lj)/u, where Jq, Ji are 



/o 1 + 

For K — oo we have Jq{uj) = 0, Jq{uj) = Jo{uj), Ji{uj) 
Bessel functions, and eq. (^ becomes the Lieb-Wu result for the ground state energy [jl3[. 
For K ^ one may rescale uj = z/k, which allows to perform the calculation analytically. 
The final resuh is E/L{U <W) = {-W + U)/A~Uy{12W), E / L{U >W) = -iyV(12f/), 
with the band- width W = 27r, in complete agreement with Ref. 0. E/L is monotonically 
increasing as function oi U/t and shows no singular behavior as function of k. 
The pseudo particle dispersions become 
du 



^'^^^ Jo cosh(c^f//4) 
e,{k) = e{k)+2 



cos{iu X)Ji{uj) + sm{ujX)J^{iu) 

'cos [iuH{k)] J^(cu) + sin [iuH{k)] J^{uj] 



(10) 



/o 1 + 

The gap A/i = /i(n 1+)— yu(n 1^) = U—2fi{n — > 1^) is given by A/i = t/— 2Max(ec(A;)). 
For large U/W one finds A/i(t/ > VT) = U-W + 0{W^/U). The analytical structure of A/x 
is very similar for all k > such that the physics cannot be different from k = oo where 
Uc = 0+. Consequently, Uc{k > 0) = 0"*", and the gap is exponentially small for U <C W. The 



situation changes for k = 0. Exphcitly, e^k, k = 0) = 2k + U - + _ 2WUk/TT /4 
such that A/i = for f/ < W, and Afi = U - W for U >W, i.e., Uc{k = 0) = W, in 
agreement with Ref. |0 and the discussion above. 

The ground state energy for strong coupling and half-filling can be cast into the form 
(2sinh(fi;)74 = W/2; iIj{x) is the digamma function) 
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E IL = / — eix) [ dx'p(x)Iie 

' U J -A 271 ^ ' J-A ' ^ ' 



e[x) = 2smh(/t) ?7 (e^ (2xsinh(K))j — r] {e_ (2xsinh(K)) 



1 X — x' 
- + i 

2 2 

-I, 



,X — X 



p{x) 



sinh(K) 



TT 



_^{2x sinh(K)) 



e' e 



(11a) 
(lib) 

(11c) 



e' (^e_^{2x sinh(/«))^ 

Here, k± = e^^{y) are the two solutions of e{k) = y with e'(A;+) > 0, e'{k ) < 0. Numerical re- 
sults for K = oo; 1; 0.1; are limu->oo{UE/L) = -4 In 2 -2.773; -2.826; -3.197; ^ 
—3.290, respectively. The expressions ([TT|) agree with the result obtained by Sutherland 
et al. [|l^ which supports our assumptions about integrability of the model ([^). They also 
provide an explicit solution for the integral equations in Ref. |1^ for the Inozemtsev model. 
We have calculated the critical exponents that control the asymptotic behavior of correlation 
functions by finite size scaling in conformal theory |T^. The final formulae generalize those 



of Ref. |T5[ to the case of different velocities for right- and left-moving electrons [|T^ . 

In sum, we are confident that the model (Q) is integrable and can thus be solved using the 
Asymptotic Bethe Ansatz. A complete construction of scattering states at finite densities 
is still lacking even for the 1/ sinh^(Kr)-Heisenberg model |P. Models with variable-range 
exchange thus remain an interesting problem in mathematical physics. 

We thank P. Nozieres and A. Gogolin for stimulating discussions. F. G. thanks all his 
colleagues at the ILL Grenoble for their kind hospitality during a stay in 1993/94. 
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